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Abstract 

We consider two strongly correlated two-component quantum systems, consisting 
of quantum mobile particles and classical immobile particles. The both systems are 
described by Falicov-Kimball-like Hamiltonians on a square lattice, extended by 
direct short-range interactions between the immobile particles. In the first system 
the mobile particles are spinless fermions while in the second one they are hardcore 
bosons. We construct rigorously ground-state phase diagrams of the both systems 
in the strong-coupling regime and at half-filling. Two main conclusions are drawn. 
Firstly, short-range interactions in quantum gases are sufficient for the appearance 
of charge stripe-ordered phases. By varying the intensity of a direct nearest-neighbor 
interaction between the immobile particles, the both systems can be driven from a 
phase-separated state (the segregated phase) to a crystalline state (the chessboard 
phase) and these transitions occur necessarily via charge-stripe phases: via a diago- 
nal striped phase in the case of fermions and via vertical (horizontal) striped phases 
in the case of hardcore bosons. Secondly, the phase diagrams of the two systems 
(mobile fermions or mobile hardcore bosons) are definitely different. However, if the 
strongest effective interaction in the fermionic case gets frustrated gently, then the 
phase diagram becomes similar to that of the bosonic case. 
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1 Introduction 



In the passing decade, some specific phases, having quasi-one-dimensional 
structure, the so called striped phases, have been a highly debated subject in 
condensed-matter physics. Apparently, a broad interest in such phases was ini- 
tiated by reports presenting experimental evidence for the existence of charge 
strips in doped layered perovskites, some of which constitute materials ex- 
hibiting high-temperature superconductivity [1,2]. 

However, striped phases had been observed much earlier, for instance on 
metallic surfaces covered with an adsorbate, see [3] and references therein. 
Theoretical descriptions of these phenomena involved various kinds of lattice 
gas models [4]. Many other instances of experimental observations of stripe- 
ordered phases are listed in [5], where Monte-Carlo studies of formation of 
striped phases in a continuous gas with hardcore and short-range repulsive 
interactions are reported. 

Quite interestingly, theoretical studies of striped phases in systems of strongly 
correlated electrons have preceded experimental observations of such phases 
[6,7]. But only after those observations, the discussion became much more vig- 
orous, and the nature of stripe-ordered phases started to attract attention of 
numerous researchers. In the context of the Hubbard model, a comprehensive 
review of the problem can be found in [8]. The existence of the same kind 
of striped phases was investigated also in the t-J model [9,10]. A bird's eye 
view on the problem of stripe-ordered phases in high-temperature supercon- 
ductors, but emphasizing its general relevance for contemporary condensed- 
matter physics can be found in [11]. The general relevance of striped phases 
is underlined also in [12], where they are viewed as an emergent phenomenon 
resulting from collective motions of microscopic particles, somewhat analogous 
to quasi-particles like phonons. Striped phases have been also studied in the 
framework of quantum-spin models, like XY model for instance [13]. 

The Hubbard or t-J-like models belong to the most realistic models, in the 
framework of which the problem of striped-ordered phases in doped layered 
perovskites can be investigated. In the both models, the spin and the charge 
degrees of freedom are taken into account, and it is believed that it is the 
competition between these degrees of freedom that is decisive for the forma- 
tion of striped phases. The question of formation of striped phases is closely 
related to the question of relative stability of these phases against mixtures 
of an electron-reach phase and a hole-reach phase. Due to the tiny energy dif- 
ferences between both phases, the results obtained by means of approximate 
methods, which introduce hardly-controllable errors, are disparate. And this 
can be said despite, as emphasized in [8], a spectacular consensus of results 
concerning half-filled stripes, obtained by means of various approximate meth- 
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ods. Therefore, as pointed out in [12,15] further careful studies are necessary 
to settle the problem of formation of stripe-ordered phases. One of possible 
ways of attacking this problem is to formulate analog problems in less realistic 
but simpler models, where some control over the results obtained by means of 
various methods of many-body physics can be gained. This leads hopefully to 
a deeper insight into the, mentioned above, stability problem. 

Such an approach has been adopted by many researchers, see for instance 
[12,14,15,16,17]. Buhler et al [14] have found, by means of Monte Carlo sim- 
ulations, that upon hole doping antiferromagnetic spin domains and charge 
stripes, whose properties are in very good agreement with experiments, ap- 
pear in a spin-fermion model for cuprates. Using the so called restricted phase 
diagrams, the stability problem of charge-stripe phases has been studied in 
the spinless Falicov-Kimball model by Lemanski et al [15,16]. In their study 
the formation of charge stripe-ordered phases can be looked upon as a way 
the system interpolates between a periodic charge- density wave phase (the 
chessboard phase) and the segregated phase (a mixture of completely filled 
and completely empty phases), as the degree of doping varies. A considerable 
reduction of the Hilbert space dimension in a spinless fermion model with in- 
finite nearest-neighbor repulsion (as compared to a Hubbard model) has been 
exploited by Zhang and Henley [12,17] to study carefully, by means of an exact 
diagonalization technique, the formation of charge stripe-ordered phases upon 
doping. They addressed also an interesting question of the role of quantum 
statistics in the problem of striped phases, by replacing fermion particles with 
hardcore boson ones. 

Our work has been inspired mainly by the recent studies of Lemanski et al 
[15,16] and by Zhang and Henley [12,17]. To investigate the key question, 
whether charge stripe-ordered phases are stable compared to a phase-separated 
state, we study rigorously the strong-coupling limit of an extended spinless 
Falicov-Kimball model on a square lattice, with mobile particles being spinless 
fermions or hardcore bosons. The usual spinless Falicov-Kimball Hamiltonian 
(such as that studied in [18]) has been augmented by a direct, Ising-like, 
interaction between the immobile particles. 

To explain the role played by the extra interaction we have to make a few 
remarks. Firstly, for technical reasons, our analysis is restricted to the case of 
half-filling. Secondly, it is now well established (proven) [18,19] that, in the 
fermion spinless Falicov-Kimball model the phase-separated state, the segre- 
gated phase, is stable only off the half-filling. Thus, an additional interaction is 
needed to stabilize the segregated phase at half-filling. Thirdly, in the regime 
of singly occupied sites and for strong-coupling, the strongest effective inter- 
action in the Hubbard model is the Heisenberg antiferromagnetic interaction. 
Doping of holes can be seen as means to weaken the strong tendency to- 
ward antiferromagnetic ordering and to make possible a phase separated state 
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with hole-reach and electron-reach regions. Under similar conditions, in the 
spinless Falicov-Kimball model an analogous role is played by the Ising-like 
nearest-neighbor (n.n.) repulsive interaction between the immobile particles 
that favors a chessboard-like ordering. The effect of weakening of the strong 
tendency toward chessboard ordering can be achieved by an extra Ising-like 
n.n. interaction that compensates the strongest repulsive interaction, and con- 
sequently permits the system to reach a phase-separated state, the segregated 
phase. On varying, in a suitable interval of values, the corresponding inter- 
action constant, which is our control parameter, we can study how the both 
systems "evolve" from the crystalline chessboard phase to the segregated one. 

We prove that by varying the interaction constant of the direct n.n. interac- 
tion between the immobile particles, the both systems (fermion and boson) 
can be driven, via a sequence of phase transitions, from a crystalline state 
(the chessboard phase) to a phase-separated state (the segregated phase). 
Moreover, these transitions occur necessarily via charge-stripe phases: via a 
diagonal striped phase in the case of fermions and via vertical (horizontal) 
striped phases in the case of hardcore bosons. Thus, short-range interactions 
in quantum gases are sufficient for the appearance of charge stripe-ordered 
phases. 

In our studies, we include also a subsidiary direct interaction between the im- 
mobile particles, an Ising-like next nearest-neighbor (n.n.n.) interaction, much 
weaker than the n.n. interaction. This interaction can reinforce or frustrate the 
n.n. interaction, depending on the sign of its interaction constant. It appears 
that on varying the control parameter (the strength of n.n. interactions), the 
systems (fermion or boson) characterized by different values of n.n.n. interac- 
tions, may undergo different sequences of transitions between the chessboard 
phase and the segregated one. We find, in particular, that if we set appropri- 
ate values of the n.n.n. interaction in the fermion system and in the boson 
one, then the both systems settle in the same phases, for typical values of 
the control parameter. On the other hand, the n.n.n. interaction enables us 
to demonstrate that the ground-state phase diagrams in the cases of mobile 
fermions and mobile hardcore bosons are definitely different. 

As a byproduct of our considerations, we obtain an Ising-like short-range 
Hamiltonian, with pair interactions at distance of two lattice constants and 
four-site (plaquette) interactions, whose ground states consist of stripe-like 
configurations or chessboard-like configurations, depending on the sign of the 
unique coupling constant. These ground states are stable with respect to per- 
turbations by n.n. and n.n.n. interactions. 

The paper is organized as follows. In the next section we present the consid- 
ered models and general properties of the corresponding Hamiltonians. Then, 
in section 3 we construct phase diagrams due to truncated effective interac- 
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tions. After that, in section 4, we analyze the phase diagrams due to complete 
interactions (without truncation) of the both quantum systems studied here 
and formulate our conclusions. Finally, in section 5, we summarize our results. 



2 The model and its basic properties 



The model to be studied is a simplified version of the one band, spin 1/2 
Hubbard model, known as the static approximation (one sort of electrons hops 
while the other sort is immobile), augmented by a direct Ising-like interaction 
V between the immobile particles. Thus the total Hamiltonian of the system 
reads: 



H 



FK 



H Q = H FK + V, 

-t 2 (4c y + C+C x ) + (c+C x - -) S x , 



T/ _ W ST ST 

v ~ g SxS y s x s y 

{x,y)i (x,y)2 



(1) 

(2) 
(3) 



In the above formulae, the underlying lattice is a square lattice, denoted A, 
consisting of sites x, y, . . ., whose number is |A|, having the shape of a \f\A] x 

\f\A~\ torus. In (2,3) and below, the sums J2{ x ,y) t > i — 1, 2, 3, stand for the 
summation over all the i-th. order n.n. pairs of lattice sites in A, with each 
pair counted once. 

The subsystem of mobile spinless particles is described in terms of creation and 
annihilation operators of an electron at site x: c+, c x , respectively, satisfying 
the canonical anticommutation relations (spinless electrons) or the commuta- 
tion relations of spin 1/2 operators S^,S~,S^, (hardcore bosons). The total 
electron-number (boson-number) operator is N e = X^c+c^, and (with a little 
abuse of notation) the corresponding electron (boson) density is p e = N e /\A\. 
There is no direct interaction between mobile particles. Their energy is due 
to hopping, with t being the n.n. hopping intensity, and due to the interac- 
tion with the localized particles, whose strength is controlled by the coupling 
constant U . 

The subsystem of localized particles (here-after called ions), is described by 
a collection of pseudo-spins {s x } xeA , with s x = 1, —1 (s x = 1 if the site x is 
occupied by an ion and s x — — 1 if it is empty), called the ion configurations. 
The total number of ions is iVj = ^2 x (s x + l)/2 and the ion density is p« = 
7Vj/|A|. In our model the ions interact directly, with energy given by V. 
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Clearly, in the composite system, whose Hamiltonian is given by (1), with 
arbitrary electron-ion (boson-ion) coupling U, the particle-number operators 
N e , JVj, and pseudo-spins s x , are conserved. Therefore the description of the 
classical subsystem in terms of the ion configurations S = {s x } xeA remains 
valid. Whenever periodic configurations of pseudo-spins are considered, it is 
assumed that A is sufficiently large, so that it accommodates an integer number 
of elementary cells. 

Nowadays, Hfk is widely known as the Hamiltonian of the spinless Falicov- 
Kimball model, a simplified version of the Hamiltonian put forward in [20]. 
To the best of our knowledge, the spinless-fermion Falicov-Kimball model is 
the unique system of interacting fermions, for which the existence of a long- 
range order (the chessboard phase) [21,22] and of a phase separation (the 
segregated phase) [18,19] have been proved. A review of other rigorous results 
and an extensive list of relevant references can be found in [23,24]). 

In what follows, we shall study the ground-state phase diagram of the system 
defined by (1) in the grand-canonical ensemble. That is, let 

H (/i e , m) = H - /j, e N e - fiiNi, (4) 

where He, Hi are the chemical potentials of the electrons (bosons) and ions, 
respectively, and let E s (/i e , Hi) be the ground-state energy of H (/i e , Hi), for a 
given configuration S of the ions. Then, the ground-state energy of H (fi e , Hi), 
Eq (He, Hi), is defined as Eq (fie, Hi) = min {Es (He, fa) '■ S}. The minimum 
is attained at the set G of the ground-state configurations of ions. We shall 
determine the subsets of the (n e , A*i)-plane, where G consists of periodic con- 
figurations of ions, uniformly in the size of the underlying square lattice. 

In studies of grand-canonical phase diagrams an important role is played by 
unitary transformations (hole-particle transformations) that exchange parti- 
cles and holes, c+c x — > 1 — c x c x , s x — » — s x , and for some (h®, Hi) leave the 
Hamiltonian H (// e , Hi) invariant. For the mobile particles such a role is played 
by the transformations: c x — > e x c x , where e x — 1 for bosons while for fermions 
e x — 1 at the even sublattice of A and e x = —1 at the odd one. Clearly, 
since H is invariant under the joint hole-particle transformation of mobile 
and localized particles, H (/i e , Hi) is hole-particle invariant at the point (0, 0). 
At the hole-particle symmetry point, the system under consideration has very 
special properties, which simplify studies of its phase diagram [22] . Moreover, 
by means of the defined above hole-particle transformations one can deter- 
mine a number of symmetries of the grand-canonical phase diagram [25] . The 
peculiarity of the model is that the case of attraction (U < 0) and the case of 
repulsion (U > 0) are related by a unitary transformation (the hole-particle 
transformation for ions): if S — {s x } xeA is a ground-state configuration at 
(He, Hi) f° r U > 0, then —S = {— s x } x£A is the ground-state configuration at 
(He, —^i) for U < 0. Consequently, without any loss of generality one can fix 
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the sign of the coupling constant U. Moreover (with the sign of U fixed), there 
is an inversion symmetry of the grand-canonical phase diagram, that is, if S 
is a ground-state configuration at (fa, fa), then —S is the ground-state config- 
uration at (— /i e , —fa)- Therefore, it is enough to determine the phase diagram 
in a half-plane specified by fixing the sign of one of the chemical potentials. 

Our aim in this paper is to investigate the ground-state phase diagrams of our 
systems, for general values of the energy parameters that appear in H (fa, fa). 
According to the state of art, this is feasible only in the strong-coupling regime, 
i.e. when \t/U\ is sufficiently small. Therefore, from now on we shall consider 
exclusively the case of a large positive coupling U, and we express all the 
parameters of H (fa,, fa) in the units of U, preserving the previous notation. 

In the strong-coupling regime, the ground-state energy E s (fa, fa) can be ex- 
panded in a power series in t. One of the ways to achieve this, for fermions 
and for hardcore bosons, is a closed-loop expansion [26,27]. The result, with 
the expansion terms up to order four shown explicitly, reads: 



E{~ (fa, fa) = -± (fa -fa)Y,s x - 1 - (fa + fa + l) | A| + 



t. - — — 
I ~ 16 + ~~8 



E s % s y + 



16 16 



E S xSy + 



t 4 



E ^ + ^E( 1 + 5^p) + j r (4) , 



{x,y)3 



(5) 



for fermions, and 



E 



boson 



(fa, fa) = ~ (fa - fa) E s x - ^ + fa + l) I A| + 



4 



X 

5^ W_ 



E s x s y 



5t 4 



e 



16 16 



E s '£ s y + 



| E ^-^E(5 + *p) + ^ (4) , 

(x,y)3 P 



(6) 



for hardcore bosons, up to a term independent of the ion configuration and 
the chemical potentials. In (5) and (6), P denotes the (2 x 2)-plaquette of 
the square lattice A, sp stands for the product of pseudo-spins assigned to 
the corners of P, and the remainders PS A \ w 4 \ which are independent of the 
chemical potentials and the parameters W and e, collect those terms of the 
expansion that are proportional to t 2m , with m — 3,4, .... It can be proved 
that the above expansions are absolutely convergent, uniformly in A, provided 
that t < 1/16 and \fa\ < 1 — I6t [26,27]. Moreover, under these conditions the 
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particle densities satisfy the half-filling relation: p e + pi — 1. 



3 Phase diagrams according to truncated expansions 

We note that on taking into account the inversion symmetry of the phase dia- 
gram in the (p e , //j)-plane and the fact that the ground-state energies depend 
only on the difference of the chemical potentials, in order to determine the 
phase diagram in the stripe \p e \ < 1 — 16t it is enough to consider the phase 
diagram at the half-line p e — 0, pi < (or p { > 0). At this half-line we set 
Pi = p. 

Due to the convergence of the expansions (5) and (6), it is possible to es- 
tablish rigorously a part of the phase diagram (that is the ground-state con- 
figurations of ions are determined everywhere in the (p e , //j)-plane, except 
some small regions) , by determining the phase diagram of the expansion trun- 
cated at the order k, that is according to the k-th order effective Hamiltonians 
^ermionyk) ( 0) ^ and ^boson^k) ( 0) ^ In order to construct a p h ase diagram 

according to a k-th order effective Hamiltonian, we use the m-potential method 
introduced in [28], with technical developments described in [25,29,27]. 

3. 1 The zeroth and second order ground-state phase diagrams 

An inspection of (5) and (6) reveals that up to the second order the effec- 
tive Hamiltonians for fermions and hardcore bosons are the same. Hence, the 
discussion that follows applies to the both cases, and the common effective 
Hamiltonians are denoted as eP (0,p), eP (0,p). 

In the zeroth order, 

E P (°,^) = ~f E (S X + 1) + y E S x S y -^ ]T S xSy 

x (x,y)i {x,y)2 

=E^ 0) , (7) 

p 

where 

h p ] = -| E' + !) + ^ E' s * s v - s ^> ( 8 ) 

and the primed sums in (8) are restricted to a plaquette P. For any e > 0, 
the plaquette potentials ifp ^ are minimized by the restrictions to a plaquette 
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w/t 1 



Fig. 1. Ground-state phase diagram of 
(0, fi) for fermion and hard-core 
boson systems. 



Fig. 2. Ground-state phase diagram of 

(2) 

E s (0, n) for fermion and hard-core 
boson systems. 



P of a few periodic configurations of ions on A. Namely, S*_ — the empty 
configuration (S + — the completely filled configuration), where s x = —1 at 
every site (s x = +1 at every site), and the chessboard configurations S^ b , 
where s x = e x , and S° b (where s x = —e x ), with e x — 1 if x belongs to the even 
sublattice of A and e x = — 1 otherwise. Moreover, out of the restrictions of 
those configurations to a plaquette P, S-\p, S + \p, S^ P , and S° b \ P , only four 
ground-state configurations on A can be built, which coincide with the four 
configurations named above. Clearly, this is due to the direct n.n.n. attractive 
interaction between the ions. The section of the phase diagram in the (W, /i, e) 
space by a plane with e > 0, according to the effective Hamiltonian (0, fi), 
is shown in Fig. 1. 

In the second order, 



4 2) (o,d 



|E(*. + i) + 



t 2 W 
4 + Y 



^2 SxS y ^2 



s x s 



{x,y)2 



(9) 



where H 



(2) 
P ) 



H 



(2) 



e_ w 

~8 + 16 



S X Sy 



" 16/ \ 

{x,y)2 



(10) 



(2) 

and the primed sums are restricted to a plaquette P. Since E s differs from 



E^P by the strength of n.n. interactions, the only effect on the phase diagram is 
that the coexistence lines are translated by the vector (— 2t 2 , 0, 0). The section 
of the phase diagram in the (W,/i,e) space by a plane with e > 0, according 
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to the effective Hamiltonian E$ ' (0,/i), is shown in Fig. 2. 



3.2 The fourth order ground-state phase diagrams 



It follows from the investigation of the phase diagram up to the second or- 
der that the degeneracy is finite, independently of the size of A, everywhere 
except the coexistence lines in the (e = 0)-plane, where it grows exponen- 
tially with |A|. Only there the effect of the fourth order interactions can be 
most significant. The meeting point of the coexistence lines in the (e = 0)- 
plane, where the coexistence line of S + and «S_ sticks to the stability domain 
of chessboard configurations, appears to be particularly interesting. In that 
point, the energies of all the configurations are the same. In what follows, 
we shall study the phase diagrams of spinless fermions and spinless hardcore 
bosons up to the fourth order, in a neighborhood of radius 0(t 4 ) of the point 
(— 2t 2 ,0,0). In this neighborhood, it is convenient to introduce new coordi- 
nates, (W,/i,e) — > (u,S,e), 

W = -2t 2 + t 4 u, n = t 4 5, e = t 4 e, (11) 
and a new (equivalent to ( y Eg ermion )^) ^-independent effective Hamiltonian, 

(E f s ermion )W (o, 5) = ^(#/™«)(4). ( 12 ) 

Then, in the spirit of the m-potential method, we express (_/j^ rTmon )(4) ^ e 
potentials ( H £ ermion )W, 

(H^ rmion )W = J2(Hl ermion Y 4 \ (13) 

T 

where, in terms of the new variables, 



{H Urmionp = _^ ^ + 1) + ^ L - ^ £" S x S y + 1 (3 - s) £" S x S y + 

{x,y)i "" (x,j/> 2 

{x,y)3 p 

In (13) and (14), T stands for a (3 x 3)-plaquette (later on called the T- 
plaquette) of a square lattice, whose sites are labeled from the left to the 
right, starting at the bottom left corner and ending in the upper right one. 
Consequently, s 5 is the pseudo-spin of the central site. The double-primed 
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Fig. 3. All the T-plaquette configurations, up to the symmetries of Hq. The config- 
urations that can be obtained from the displayed ones by the hole-particle transfor- 
mation are not shown, only the numbers assigned to them are given in brackets. 

sums are restricted to a T-plaquette. For bosons, we introduce (H^ son )^ and 
(H^ son )W in the same manner as for fermions, with 



(i4° s ° n ) (4) = S (*b + 1) + (w - I) E" ***v + ^ (5 - e) E" + 

-^Y!' s * s y- ^E"(^ + 5 )- ( 15 ) 



By the remark at the begining of this section, we have to search for the lowest- 
energy configurations of (^/ errmon )(4) anc i ^boson^(4) amorL g a }[ ^he config- 
urations. Consequently, the potentials ^f/// erm * 0?1 )(4) anc i ^jjboson-j{4) have to 
be minimized over all the T-plaquette configurations. There are (up to the 
symmetries of H ) 102 different T-plaquette configurations, shown in Fig. 3. 
Unfortunately, in contrast to the lower-order cases, the potentials (f/"/ errmon )(4) 
and (if|? son )( 4 ) turn out to be the m-potentials only in a small part of (u, 5, e)- 
space. This difficulty can be overcome by introducing the so called zero- 
potentials [25,29,27], denoted that are invariant with respect to the 
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symmetries of H and satisfy the condition 



E4 4) = 0. (16) 



Following [29], are chosen in the form 



4 4) = ^>4 i} , (17) 
i=i 



where the coefficients ai, depending on (u,5,e) in general, have to be deter- 
mined in the process of constructing the phase diagram, and the potentials 
left are defined by 



k^ = si + s 3 + s 7 + s 9 - 4s 5 , 

k y T — S 2 + S4 + Sq + Sg - 4s 5 , 
(3) 

k V T ' = SiS 2 + S 2 S 3 + S 3 S 6 + S e Sg + S 8 S 9 + S 7 Sg + S 4 S 7 + S1S4 



— 2s 2 s 5 — 2S5S6 — 2s 5 s 8 — 2S4S5, 
k^ = S1S3 + S3S9 + S7S9 + sis 7 — 2s 4 s 6 — 2s 2 s 8 . (18) 



, (4) 

k K T ' = S1S5 + s 3 s 5 + S5S9 + s 5 s 7 - S2S4 - s 4 s 8 - s 8 s 6 - s 2 s 6 , 



Note that the potentials k^ are invariant with respect to the symmetries of Hq 
and they satisfy (16), therefore these properties are shared by the potentials 

4 4) . 

By means of the zero-potentials \ new candidates for m-potentials can be 
introduced, 

= + K ^ ; ( 19 ) 

T 

with an analogous representation of {H^ son )^ A \ and now the task is to min- 
imize the potentials ^ ermion )^ + K^and (H b T oson )^ + over all the 
T-plaquette configurations. 

In our study of the ground-state phase diagrams we limit ourselves to the (5 = 
0)-plane, where the both considered systems are hole-particle invariant, and to 
the (e = 0)-plane. Our analysis of the minima of the T-plaquette potentials, 
^/ er-rmon )(4) and (iJ^ oso ™)( 4 ) augmented by the zero-potentials shows 
that these planes are partitioned into a finite number of open domains Sd, 
each with its unique set of periodic ground-state configurations on A, denoted 
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also by Sd- There is a finite number, independent of |A|, of configurations in 
Sd and they are related by the symmetries of H . The last two statements do 
not apply to only one of the domains, denoted Sd2, which will be described in 
the sequel. 

The domains Sd are characterized as follows: at each point p (p = (ou,e) or 
p = (uj,5)) of a domain Sd, there exist a set of coefficients {ai(p)} such that 
the corresponding potentials are minimized by a set Std(p) of T-plaquette 
configurations. Moreover, from the configurations in Std(p) one can construct 
only the configurations in S D - The set of the restrictions to T-plaquettes of 
the configurations from Sd, Sd\t, is contained in each set Std(p) with p e Sd- 
In Table 4 of Appendix we mark by asterisk the cases, where the set Std{p) 
contains, besides Sd\t, some additional T-plaquette configurations. 

Specifically, for 5 = the ground-state phase diagrams due to the effective 
Hamiltonians (H^ rmion )^ and (H^ son )W are shown in Fig. 4 and Fig. 5, 
respectively, while the corresponding ground-state phase diagrams for 5 = 0, 
in Fig. 6 and Fig. 7. In Fig. 8 we display the representatives of the sets Sd 
of ground-state configurations. That is, the remaining configurations of Sd 
can be obtained easily by applying the symmetries of H to the displayed 
configurations. The domain S seg , having no representative in Fig. 8, consists 
of the two translation- invariant configurations S+ and S-, related by the hole- 
particle transformation. 

Only in the domain Sd2, which appears in the phase diagrams shown in 
Fig. 6, 7, that is off the hole-particle symmetry plain, the situation is not 
that simple. In Sd2 one can distinguish two classes, Sd2a and Sd2b, of periodic 
configurations with parallelogram elementary cells. A configuration in Sd2a 
consists of vertical (horizontal) dimers of filled sites that form a square lat- 
tice, where the sides of the elementary squares have the length 2\[2 and the 
slope ±1. In a configuration of Sd2b, the elementary parallelograms formed by 
dimers have the sides of the length 2 a/2 and the slope ±1, and the sides of 
the length vTO and the slope ±1/3. Two configurations, one from Sd2a and 
one from Sd2b, having the same kind of dimers (vertical or horizontal), can be 
merged together along a "defect line" of the slope ±1 (dashed line in Fig. 8), 
as shown in Fig. 8, without increasing the energy. By introducing more defect 
lines one can construct many ground-state configurations whose number scales 
with the size of the lattice as exp {const \fh). While for a finite lattice all the 
configurations are periodic, many of them become aperiodic in the infinite 
volume limit. 

As mentioned above, except Sd2 all the other sets Sd contain exclusively pe- 
riodic configurations. The set S pc i> (of plaquette-chessboard configurations) 
contains configurations built out of elementary plaquettes with filled sites, 
forming a square lattice according to the same rules as filled sites form a 
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56,57,63,64,69,75,82,83,92,94,100} 57,63,67,68,72,78,82,84,89,94,95} 



Fig. 4. The (/x = 0)-phase diagram of the effective Hamiltonian ^jjf^ rrnion ^) m xhe 
numbers in curly brackets, displayed by the symbols of open domains So, denote 
the T-plaquette configurations that minimize the m-potential in Sd- The num- 
bers in curly brackets, displayed by the boundary-line segments or by the arrows 
pointing towards boundary segments (or their crossing points) identify the addi- 
tional minimizing T-plaquette configurations. For more comments see the text. The 
boundary-line segments can be determined by means of their intersection points: 
a = (9/2,12), b = (-5/2,-2), c = (23/2,-2), d = (-5/2,-6), e = (23/2,-6), 
and by the slope 1 of the boundary of S v i, and the slope —1 of the boundary of 5^1- 




Fig. 5. The (ji = 0)-phase diagram of the effective Hamiltonian (#&° so ")( 4 ) . The 
boundary-line segments can be determined by means of their intersection points: 
a = (5/2,9), b = (1/2,5), c = (9/2,5), d = (1/2,1), e = (9/2,1), and by the 
slope 1 of the boundary of S v i, and the slope —1 of the boundary of Sdi- For more 
explanations see the description of Fig. 4. 
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Fig. 6. The (e = 0)-phase diagram of the effective Hamiltonian (jy£ rmion )( 4 ). The 
boundary-line segments can be determined by means of their intersection points: 
a = (1/2,-1), b = (7/2,-3/2), c = (15/2,-3/2), and by the slope -1 of the 
boundaries of Si, S2, S3. For more comments see the description of Fig. 4. 




Fig. 7. The (e = 0)-phase diagram of the effective Hamiltonian (i^g so ™)( 4 ). 
The boundary-line segments can be determined by means of their intersection 
points: a = (-5/2,-1/2), b = (11/2,-1), c = (7/2,-2), d = (3/2,-5/2), 
e = (15/2,-5/2), f = (19/2,-3), g = (11/2,-7/2), h = (7/2,-4), and by the 
slope —1 of the boundaries of Si and S3. For more comments see the description of 
Fig. 4 
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5 v5 r P =7/5; S 2 (p=l/4) S d2 (p=l/4) SjfprVS) 

Fig. 8. Representative configurations of the sets Sd of ground-state configurations. 
The remaining configurations of Sd can be obtained by applying the symmetries 
of Hq to the displayed configurations. As a representative configuration of 5^2, we 
show a configuration with one defect line (the dashed line). The continuous line is 
a guide for the eye. For more comments see the text. 

square lattice in the chessboard configurations from <S C (,. The remaining sets 
Sd consist of configurations that have a quasi-one- dimensional structure. That 
is, they are built out of completely filled and completely empty lattice lines 
of given slope. Such a configuration can be specified by the slope of the filled 
lattice lines of the representative configuration and the succession of filled (f) 
and empty (e) consecutive lattice lines in a period. For instance, the repre- 
sentative configuration of Sdi (see Fig. 8) is built out of filled lines with the 
slope 2 and, in the period, two consecutive filled lines are followed by two 
consecutive empty lines, which is denoted (2;2/, 2e). Similar description of 
the remaining quasi-one- dimensional configurations is given in Fig. 8. 



The coefficients {a,} for which the fourth order potentials, (H? 



fermion\ (4) 



+ K. 



(4) 
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and (H^ oson )^ + Kf \ become m-potentials are given in the tables collected 
in the Appendix. 

Finally, a remark concerning ground-state configurations at the boundaries 
between the open domains Sd is in order. Let Sd and Sd> be two domains 
of the considered phase diagram, sharing a boundary. At this boundary, the 
set of the minimizing T-plaquette configurations contains always the subset 
Std^Std 1 , but it may contain also some additional T-plaquette configurations 
of minimal energy. Consequently, the set of the ground-state configurations at 
the boundary contains always the subset Sd USd' and typically a great many 
of other ground-state configurations, whose number grows indefinitely with 
the size of the lattice. In the considered diagrams, the only exception is the 
boundary between S seg and S c b, where the set of the ground-state configura- 
tions amounts exactly to S seg U S c b- 



4 Discussion of the phase diagrams and conclusions 

We start with a few comments concerning the phase diagrams due to the 
fourth order effective interactions. In the fourth order effective interaction, 
the parameters uj and e control the strength of n.n. and n.n.n. interactions, 
respectively. The n.n. interaction is repulsive and favors the chessboard con- 
figurations if uj > 9/2, for fermions (uj > 5/2, for bosons). In the opposite case 
(attraction) it favors the S + , S_ configurations. In turn, the n.n.n. interaction 
is attractive if e > 3, for fermions (e > 5, for bosons), and then it reinforces 
the tendency towards chessboard and uniform, S + , S_, configurations. When 
it becomes repulsive, it frustrates the n.n. interaction: the more negative it 
is the larger uj is needed to stabilize the chessboard configurations and the 
smaller uj is needed to stabilize the S + , S_ configurations. Therefore, what- 
ever the value of e is, there is a sufficiently large uj such that (e, uj) G S c b and 
a sufficiently small uj such that (e,u) G S seg . 

Apparently, the fourth order phase diagrams in the cases of hoping fermions 
and hoping bosons are quite similar if we compare the geometry of the phase 
boundaries and the sets of ground states. The main difference is in the domain 
occupying the central position: in the case of fermions the ground-state config- 
urations are diagonal-stripe configurations, Sdd, while in the case of bosons this 
is the set S pc b of plaquette-chessboard configurations. These two sets of con- 
figurations are ground-state configurations of the corresponding fourth-order 
effective interactions with the site, n.n., and n.n.n. terms dropped, which cor- 
responds to the points (9/2, 3) G Sdd and (5/2, 5) G S pc b- This means that the 
fourth order phase diagrams can be looked upon as the result of perturbing 
an Ising-like Hamiltonian, that consists of pair interactions at distance of two 
lattice constants and plaquette interactions, by n.n. and n.n.n. interactions. 
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More generally, it is easy to verify that the set of ground-state configurations 
of the Ising-like Hamiltonian, 

J2 SxSy + "fY. S P> ( 20 ) 

(x,y) 3 P 

amount to S<id if 7 > 0, and to S pc i, if 7 < 0. These ground states are stable with 
respect to perturbations by n.n. and n.n.n. interactions, if the corresponding 
interaction constants are in a certain vicinity of zero. Otherwise, new sets of 
ground states, shown in Fig. 4 and Fig. 5 emerge. 

The basic question to be answered, before discussing the phase diagrams due to 
the complete interaction, refers to the relation between these phase diagrams 
and the diagrams due to the truncated effective interactions, obtained in the 
previous section. 

Firstly, we note, by inspection of the phase boundaries in Fig. 4 and Fig. 5, 
that the T-plaquette configurations {45, 58} are missing in the phase diagram 
corresponding to fermion mobile particles. In turn, the T-plaquette config- 
urations {41,62} are missing in the phase diagram corresponding to boson 
mobile particles. Consequently, the ground-state configurations of S pc i, can- 
not appear not only in the fourth order fermion phase diagram but also in the 
fermion phase diagram of the complete interaction. Similarly, the ground-state 
configurations of Sm do not appear in boson phase diagrams. 

Secondly, by adapting the arguments presented in [29,27], we can demonstrate, 
see for instance [31], that if the remainder is taken into account, then there 
is a (sufficiently small) constant to such that for t < t the phase diagram looks 
the same as the phase diagram due to the effective interaction truncated at the 
fourth order, except some regions of width 0(t 6 ), located along the boundaries 
between the domains, and except the domain Sd2- For t < to and each domain 
Sd, <Sd 7^ <Sd2, there is a nonempty two-dimensional open domain that 
is contained in the domain Sd and such that in the set of ground-state 
configurations coincides with S D . 

Now, consider our systems for specified particle densities, p e = Pi = 1/2. Ac- 
cording to the phase diagrams in the hole-particle symmetry plane (p = 5 = 
0), shown in Fig. 4 and Fig. 5, if (e,uj) is in S™ g , then the ground state is 
a phase-separated state, where ion configurations are mixtures of S + and 5_ 
configurations, called the segregated phase. Another phase-separated state, 
where the ion configurations are mixtures of S + , S-, S^ b , and S° b configura- 
tions, exists at a line that is a small (0(t 2 )) distortion of the boundary between 
the domains S seg and S c b [29]. In all the other domains of these diagrams 
the ground-state phase exhibits a crystalline long-range order of ions. 
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By the properties of the fourth order phase diagrams, whatever the value of e 
is, there is a sufficiently large oo such that the systems are in the chessboard 
phase ((e,u) G S°g) and a sufficiently small u such that the systems are in 
the segregated phase ((s,cu) G S™ g ). In between, the systems undergo a series 
of phase transitions and visit various striped phases. Consider first unbiased 
systems (e — 0). Then, it follows from the phase diagram in Fig. 4, that the 
fermion system visits necessarily the diagonal striped phase, described by the 
configurations in Sdd- In turn, the phase diagram in Fig. 5 implies that the 
boson system has to visit the dimerized-chessboard phase, described by the 
configurations in iS^i, then the vertical/horizontal striped phase, described by 
the configurations in S V 2, and after that a dimerized vertical/horizontal striped 
phase, described by the configurations in S v i. These scenarios are preserved, 
if e G (-2 + 0(t 2 ), 12 - 0{t 2 )) in the fermion case, and e G (-oo, 1 - 0{t 2 )) 
in the boson case. 

Note that the same sequence of transitions, as found in the unbiased boson 
system, can be realized by the fermion system if e is sufficiently small (e < 
-6-0(t 2 )). 

For somewhat larger values of e, e G (—6 + 0(t 2 ), —2 — 0(t 2 )), on its way 
from the chessboard phase to the segregated phase the fermion systems visits 
the dimerized chessboard phase, then the diagonal striped phase, and after 
that the dimerized vertical/horizontal striped phase. The vertical/horizontal 
striped phase is closer to the segregated phase than the diagonal striped phase. 
In case the both kinds of stripe phases (vertical/horizontal and diagonal) ap- 
pear in a phase diagram, such a succession of striped phases is perhaps generic, 
see also [15,32]. 

For sufficiently large e, the tendency toward the chessboard and the uniform 
configurations is so strong that the both systems jump directly from the chess- 
board phase to the segregated phase. 

Above, we have described the states of the considered systems for typical 
values of the control parameter uo. Only in small intervals (whose width is of 
the order of 0(t 2 )) of values of u, about the transition points of the diagrams 
in Fig. 4 and Fig. 5, the states of the systems remain undetermined. 

For completeness we present also ground-state phase diagrams of the both 
systems off the hole-particle symmetry plane, for e — 0, see Fig. 6, 7. Away 
from the hole-particle symmetry plane new phases appear. In particular we 
find the phases described by the quasi-one-dimensional configurations <Si, <S 2 , 
and S3, well known from the studies of the phase diagram of the spinless 
Falicov-Kimball model [27]. In the boson phase diagram, Fig. 7, there appears 
the domain Sd2 that is not amenable to the kind of arguments used in this 
paper. 
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5 Summary 



In this paper, we address the problem of the stability of charge-stripe phases 
versus the phase separated state — the segregated phase, in strongly-interacting 
systems of fermions or hardcore bosons. There are numerous works devoted to 
this problem, where it is studied, by approximate methods, in the framework 
of relevant for experiments models. Unfortunately, due to tiny energy differ- 
ences involved, it is difficult to settle this problem by means of approximate 
methods, which bias the calculated energies with hardly controllable errors of 
various nature. We have studied simple models, that by many physicists can 
be considered less realistic, but which, in return, are amenable to a rigorous 
analysis. Our models are described by extended Falicov-Kimball Hamiltoni- 
ans, with hoping particles being spinless fermions or hardcore bosons. The 
ground-state phase diagrams of these models have been constructed rigor- 
ously in the regime of strong coupling and half-filling. In the both cases, of 
fermions and hardcore bosons, we have found transitions from a crystalline 
chessboard phase to the segregated phase via striped phases. 
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6 Appendix 

In this section, at each point p = (ou,e) we give the sets {ai(p)} of the zero- 
potential coefficients o;, that we used to construct the four phase diagrams, 
presented in the previous section. Specifically, for each of the four phase di- 
agrams we define a partition, called a-partition, of the (uj, e)-plane (see the 
tables displayed below) into sets that differ, in general, from the domains Sd- 
To each set of such a partition we assign a set {(aj,&j,Cj)} of triplets of ra- 
tional numbers, such that the coefficients Oiiip) (with p in the considered set) 
are affine functions of the form a^u; + biE + q. We have not been able to find 
one set {(aj, bi, q)}, that turns the fourth order potentials into m-potentials 
in the whole (uj, e)-plane. We have not also succeeded in assigning to each 
domain Sd an exactly one set {(aj, bi, q)}. The same remarks apply to the 
(uj, <5)-plane. 
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Fig. 9. The sets A, B, C, whose boundaries are marked with dashed lines, used in 
Table 1 to define the a-partition of the fourth order phase diagram in the case of 
fermions and for 5 = O.The dashed-line segments are determined by their intersec- 
tion points: a = (-1/2, -6), b = (-13/2, -10), c = (19/2, -6), d = (31/2, -10). 
The coordinates of the remaining points can be read from Fig. 4. 




Fig. 10. The sets A, B, C, whose boundaries are marked with dashed lines, used in 
Table 2 to define the a-partition of the fourth order phase diagram in the case of 
hardcore bosons and for 5 = 0. The dashed-line segments are determined by their 
intersection points: a = (1/2,4), b = (2,1), c = (-5/2,-2), d = (9/2,4), e = (3,1), 
f = (15/2,-2). 
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Table 1 

The set of zero-potential coefficients {a{\ in the case of fermions and for 5 = 0. In 
the first column the sets of the a-partition are specified. For more comments see 
the text in Appendix. 
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Table 2 

The set of zero-potential coefficients {c^} in the case of hardcore bosons and for 5 = 
0. In the first column the sets of the a-partition are specified. For more comments 
see the text in Appendix. 
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Table 3 



The set of zero-potential coefficients in the case of fermions and for e = 0. In 
the first column the sets of the a-partition are specified. For more comments see 
the text in Appendix. 
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Fig. 11. The sets A,B,. . . , whose boundaries are marked with dashed and continuous 
lines (the continuous lines are the boundaries of the phase diagram), used in Table 3 
to define the a-partition of the fourth order phase diagram in the case of fermions 
and for e = 0. 

-6 -4 -2 2 4 6 8 10 12 
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Fig. 12. The sets A,B,. . . , whose boundaries are marked with dashed and contin- 
uous lines (the continuous lines are the boundaries of the phase diagram) , used in 
Table 4 to define the a-partition of the fourth order phase diagram in the case of 
hardcore bosons and for e = 0. The dashed-line segments are determined by their in- 
tersection points: a = (-5/2, -1), b = (3/2, -3), c = (7/2, -13/3), d = (9/2, -5), 
e = (6, -4), f = (29/4,-3). 
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Table 4 

The set of zero-potential coefficients {a{\ in the case of hardcore bosons and for 
e = 0. In the first column the sets of the a-partition are specified. The cases, where 



the set Sd\t 1S a proper subset of Std(p) are marked by the asterisk. For more 
comments see the text in Appendix. 
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20 20 ' 4 


uj 4 5 
60 20 96 


uj 4 23 
32 16 64 


3uj 34 43 
160 80 192 


«S d3 nF * 


4 

16 


3w 34 , 5 
160 20 + 64 


31uj 4 , 13 
960 20 + 384 


4 3 
16 16 


uj 34 53 
320 80 384 


«S d3 nG * 


3uj 4 1 3 
64 ' 16 " r 128 


ui 5 1 3 
16 4 + 32 


uj , 4 _i_ 5 
240 ' 20 "t" 48 


3uj _i_ 4 21 
64 ' 16 128 


9uj , 34 41 
320 ' 80 384 


«S d3 nH * 


w , 4 , 25 
64 "f" 16 ' 128 


uj 4 5 
32 4 64 


19uj , 4 , 73 
960 20 "f" 384 


uj i 4 _,_ 1 
64 16 t 128 


uj , 35 1 
80 "f" 80 48 
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